GENERAL TREATMENT OF VORTICAL, TOROIDAL, AND COMPRESSION 

MODES 
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The multipole vortical, toroidal, and compression modes are analyzed. Following the vorticity 
concept of Ravenhall and Wambach, the vortical operator is derived and related in a simple way 
to the toroidal and compression operators. The strength functions and velocity fields of the modes 
are analyzed in '^"^Pb within the random-phase-approximation using the Skyrme force SLy6. Both 
convection and magnetization nuclear currents are taken into account. It is shown that the isoscalar 
(isovector) vortical and toroidal modes are dominated by the convection (magnetization) nuclear 
current while the compression mode is fully convective. The relation between the above concept of 
the vorticity to the hydrodynamical vorticity is briefly discussed. 



I. INTRODUCTION 



An irrotational character of nuclear flow is a basic 
assumption in collective nuclear dynamics [ll-0|, which 
manifests itself in numerous examples of low-energy ex- 
citations and giant resonances (GR). At the same time, 
nuclear motion can also carry a vorticity, i.e. a deviation 
from irrotational flow In hydrodynamics (HD), 

the vorticity is defined as a curl of the velocity field [3. 
Instead, the nuclear theory deals with nuclear currents 
rather than velocities and so here the vorticity is often 
defined through the j\\^i(r) component of the multipole 
decomposition of the transverse nuclear current The 
component jx\+iif) is treated as unrestricted by the con- 
tinuity equation (containing the current divergence V ■ j) 
and so is believed to be of a vortical character. In this 
case, the vorticity and charge transition density represent 
two independent parts of the charge-current distribution. 
This consideration reminds the previous result [l3| where 
the current component jAA-i('') is proposed to be con- 
strained by the continuity equation and thus determined 
by the charge distribution while the component jxx+iif) 
is treated as independent. 

Both definitions of the vorticity, from HD and Ref. , 
are widely used in the literature: the former in the nu- 
clear fluid-dynamical models (see, e.g. [ll| and references 
therein) and the latter in the microscopic studies, see e.g. 
0, [3 . These two definitions arc assumed to be closely 
related , though they are different obscrvables by con- 
struction. Actually, they represent different aspects of 
the nuclear vorticity. In this paper, we will concentrate 
on the jAA+i ('')-based vorticity |7|. The comparison with 
the HD case will be also done. 

The most remarkable manifestation of vorticity is the 



electric dipole toroidal mode (TM) jl2l - ll4l | intensively ex- 
plored during the last decades, see e.g. the review [l5| 
and references therein. This exotic mode is treated as 
a vortical collective motion of the toroidal type. The 
TM operator is the second-order correction to the lead- 
ing El operator in the long-wave approximation. An- 
other kind of the second-order El flow is represented by 
the anisotropic compression mode (CM), often called the 
isoscalar dipole GR [T6l4l9l |. The mode is viewed as a 
compression wave in a definite direction and so is related 
to the nuclear incompressibility [l^, [l^. The isoscalar 
(T=0) TM and CM were observed in (a, a')-reaction as 
broad low-energy (TM dominated) and high-energy (CM 
dominated) electric dipole distributions [13, [H, |20[ - [25| . 
The TM was also investigated in the region of the pygmy 
resonance in ^"^^Pb in a nuclear fluorescence experiment 
[2^. Perhaps, the TM can be discriminated in the trans- 
verse (e, e') form-factors 
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The TM and CM were examined in various models, 
includin g th e fluid-dynamical and sum-rule approaches 
[H [Iljllil, the method of Wi gncr function mo- 
ments 29r 3Cl|. the random-phase-approximation (RPA) 
[3l| and multi-phonon (26j methods with phcnomcnolog- 
ical single-particle potentials. More refined RPA studies 
within the self-consistent mean-field approaches were also 
performed, relativistic ones [13, HI] and those based on 
Skyrme forces [13, HE] i for a review see [lH| • A direct 
relation between the current-dependent TM and density- 
dependent CM operators was established in [ij, |ll| and 
both modes were shown to be mixed. Most of the 
studies reproduce the observed bimodal (low-energy TM 
and high-energy CM) distribution. However, theoretical 
models generally overestimate the CM peak energy by 
--4 MeV and underestimates the TM one by 1-2 MeV. 
Besides, they yield a much broader TM distribution [isj . 

A special effort was devoted to the nuclear vorticity as 
such j^-dl- The possibility to measure the vorticity in 
(e, e') experiments was discussed [1, [27l. l36j. 

Despite these thorough studies, some principle points 
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concerning the vorticity and related modes deserve fur- 
ther inspection: i) There is an essential difference in 
modeling the vortical mode (VM) and their TM and CM 
counterparts. The TM and CM are usually deduced as 
second order terms in a low-momentum expansion of the 
electric multipole transition operators [ij, ■ To the 
best of our knowledge, an analogous way to the VM op- 
erator has yet to be developed. This would help to es- 
tablish a formal relation between VM and TM/CM. ii) 
Most of the previous studies (with exception of 0,11, |3l|) 
employ only the convection part jc of the nuclear current 
and skip its magnetization (spin) part jm, though the 
latter can also produce a vorticity. The role of jm in VM 
and other modes has yet to be clarified, iii) Mainly the 
T=0 channel of the modes were discussed, although their 
T=l counterpart is also interesting and deserves a closer 
look, iv) The relation between two definition of the vor- 
ticity, from HD Q and jxx+i{r) current component 0, 
has yet to be clarified. 

The aim of the present study is to explore the open 
problems listed above. First of all, the vortical operator 
unconstrained by the charge conservation is derived fol- 
lowing the ideas 0|- The operator has a simple relation 
with its TM/CM counterparts and is also deduced as a 
second-order term in a low-momentum (long wavelength) 
expansion of the dominant electric operator. Further, 
the difference in the vorticity criteria from the HD 
and Ref. is inspected. In the numerical calculations, 
both T=0 and T=l channels of VM, TM, and CM are 
analyzed by using the full nuclear current j = jc + jm ■ 
The dominant role of jm and thus the spin vorticity in 
forming the isovector VM and TM is worked out. Both 
single-particle and collective mechanisms of the vorticity 
are discussed. 

The numerical analysis is done within the self- 
consistent separable random-phase-approximation 
(SRPA) approach based on the factorized Skyrme 
residual interaction [s^, [11] . The systematic studies of 
electric [37l - l42| and magnetic (43l - l45| GR in spherical 
and deformed nuclei have established this approach as a 
reliable and effective theoretical tool. 

The paper is organized as follows. In Sec. II, the 
nuclear vorticity is discussed in context of the previous 
studies. In Sec. Ill, the VM, TM, and CM operators are 
derived on the same theoretical grounds, following the 
prescription 0. A simple relation between the opera- 
tors is established. The mode velocities are analyzed and 
different criteria of the vorticity, from and HD, arc 
compared. Sec. IV provides an outline of the calculation 
scheme within the Skyrme SRPA. In Sec. V, the numer- 
ical results for the T=0 and T=l VM, TM, and CM are 
discussed. A summary is given in Sec. VI. Appendix 
[^justifies the procedure of building the vortical opera- 
tor. Appendix [B] describes the procedure for extraction 
of the spurious ccnter-of-mass corrections. Appendix [C] 
provides the explicit expressions for the density and cur- 
rent operators. Appendix ID] sketches the basics of SRPA. 



II. VORTICITY IN TERMS OF VELOCITIES 
AND CURRENTS 

The HD nuclear models, including the famous liquid- 
drop model, assume an irrotational character of the col- 
lective nuclear flow [iHll 

V X v{f) = (1) 

where v{'r) is the local velocity field. The deviation from 
the irrotational flow is quantifled by the HD vorticity 

ro(r) = V X v{f) . (2) 

Unlike the HD models, nuclear theory prefers to deal 
with currents jnuc{'r) rather than velocities v{r). How- 
ever, Eqs. ((H) and © cannot be replaced by the similar 
expressions for the nuclear current since, as shown be- 
low, the curls of velocity and current have a different 
structure. Moreover, V x jnuc{r) is the key part of the 
electrical multipole operator M{EXfi, k) and, so, treating 

^ X Jnuc 

(f) as a vortical quantity would wrongly mean a 
fully vortical nature of any EXjjL excitations, which con- 
tradicts, e.g., a predominantly irrotational character of 
electric GR. 

For using the HD definition of vorticity ([5]), the quan- 
tum theory should express it through the nuclear current. 
A common way is to define the velocity field Vu {r) for the 
excitation mode v through the current transition density 
5j^{r) and ground state density po{r) [gI. [ioj. 

5jy{r) = po{r)v^{r) , (3) 
which casts ^ into the form 

V X ^ P'^"^^ ^ ^>(r) - Vpo(r) X 5U^ 

This illuminates the difference between curls of the ve- 
locity and current and thus shows that V x 5ji^ (r) cannot 
be a measure of vorticity. The difference is comprised in 
the gradient field V po{r) which is maximal at the nuclear 
surface and minimal in the interior. 

In , a concept of nuclear vorticity, alternative to the 
HD one and fully based on the nuclear current, was pro- 
posed. It aims to find a component of the nuclear current 
which is not restricted by the continuity equation 

V • 5ju{r) = -ikc5p^{r) (5) 

with k being the transfer momentum. By construction, 
this current component should not contribute to the ir- 
rotational motion and vanish in the divergence V-(5ji/(r). 
So it may be naturally used for building the vortical 
quantities. 

Since the present study follows similar lines, it is worth 
to outline the concept [7| in more detail. In spherical 
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nuclei, it exploits the multipolc expansion of the nucleon 
and current transition densities: 



SpMr) = {jfmf\p{f)\j,mi) 

Xfj, 



(6) 



E <JAL(r)fALM« , (7) 



\fi L=X±1 



where 



fi 



= (-1) 



(8) 



jii {jfy TTT-f) a-re spins and their projections for the ini- 
tial i (final /) state, Y\^i{'r) and YxL^i{r) are ordinary and 
vector spherical harmonics (47j . The p\{r) and jxL{r) 
are multipole components of the transition density and 
current. Using the above expansions and the quantity 



A + l 



kcpxir) Yxx^ir) , (9) 



the unconstrained vortical transition density 

5wfi{r) = V X S'jfiir) - 6Sf^{f) 



(10) 



is built , where the vortical multipoles 



wxxir) 



2A + 1 



A 



d 
dr 



A + 2 



r 



jxx+i{r) (11) 



are determined by the radial current component 
j\ A+i(?')- Finally, the vorticity strength is given by 



dr r'^^'* wa xii") 



(12) 



The continuity equation ((S)) in terms of the moments 
relates the current component jx x-iij") to the transition 
density but leaves the component jAA+i('') untouched 
01 • So, just jAA-i('') provides the charge-conservation 
constraint and the quantity ^ is constructed so as to 
remove jAA-i('") from V x Sjfi(r), see more discussion 
in Appendix]^ The unconstrained vortical value (jlip in- 
cludes only jA A-i-i('^)- Moreover, following Q, the motion 
is treated as vortical if its current involves jx x+iir)- 

The formalism treats the vorticity without using 
an explicit vortical operator analogous to the TM and 
CM ones. However, we need such operator for reasons of 
comparison of vortical modes with TM and CM. In the 
next section, we will develop the explicit vortical opera- 
tor by using the prescription and relate this operator 
with the toroidal and compression ones. The subsequent 
discussion will show that this operator is not truly vorti- 
cal in the HD sense and the presence of jx x+iir) in the 
current is not enough for the HD vorticity. 



III. VM, TM, AND CM OPERATORS 



A. Derivation of VM, TM and, CM operators 



The standard electrical multipole operator may be 
written in different forms 1481: 



M{EXp, k) 



(2A + 1)!! 



(13) 



cP+i(A-M) 

rfVJ„„e(r) • [V X (f X V) jA (/c^i^AM (^] 



(2A + 1)!! / A 



(14) 



cfc^+i V A + 1 
d\[Jx{kr)Yxx^.{^]■ [V x j;„,(f)] 



where jx{kr) is the spherical Bessel function. 

The form (fT4|) shows that V x jnuc{r) cannot be a mea- 
sure of the vorticity since otherwise M{EXfi, k) would 
indicate only vortical electric excitations. At the same 
time, the form <\1A\ suggests that the vortical operator 

may be built from M{EkXfi) by replacing V x jnuci^) 
with the truly vortical quantity [7| 



V X jnucir) - jkc [yp{r) x r\ 



(15) 



The density-dependent term in ([T51) subtracts the charge 
conservation constraint. Actually it plays a similar role 
as the r.h.s. second term in the vortical transition den- 
sity (|T0)) . Both prescriptions, (|10p and (|T5|) . have the 
same intention but, being applied to to different quan- 
tities (transition densities and operators), lead to for- 
mally different recipes. While ([TU]) excludes the ja a-i('') 
terms, the recipe (|15p gives an exact compensation of the 
lowest-order fc-tcrms in the final vortical operator, see the 
derivation below. Both ^TU\\ and ([T5|) remind the r.h.s. of 
the HD expression ^ . They are compared and discussed 
in more detail in Appendix [X] 

By using the vortical operator is defined as 
M^or{EXp, k) = M{EXp, k) - Ms{EXp, k) (16) 

i.e. as a difference of the electric operator and the 
subsidiary operator 



MsiEXfi,k) 



(2A+1)!! 



k^^XiX + 1) 

d\[jxikr)Yxx>.ir^]- [Vp(f) x r] 



(17) 
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The latter may be also written in the forms 



(2A + 1)!! 



(2A + 1)!! 



d^rp{r)jx{kr)Yx^{r) 



(19) 



The form (j20p is obtained by using the operator continu- 
ity equation 



is the familiar electrical operator in the long- wavelength 
limit (lowest order term) and 



(26) 



2c(A + l)(2A + 3) 



^^^^ j d'rj,,uc{r) ■ rVn^(r) (27) 
S(a/i);2A + 3) /-^^-^(->^"^-^(-^' 



V • jnuc = - - [i?, p] = -ikcp . 
n 



(21) 



In the long- wavelength approximation (k — > 0), we 
keep only the first and second terms in the expansion 
of the spherical Bessel function 



jA(fcr) 



{krY 



(kr) 



(22) 



(2A + 1)!! 2(2A + 3) ' 

and thus get for the electric and subsidiary operators 

M{EXp,k) w M{EXp) + kMtor{EXp) , (23) 
MsiEXp,k) « M{EXp)~kM,omiEXp) , (24) 

where 



MiEXfj.) = — / d'rj^^,if)-V{r^Yx^,ir)) (25) 



d^r (V-Jn«c(r))r^rA^(^ 
= - / d^rp{r)r^Yx^ir^ 



McomiEXp) 



2c(2A + 3) 



(28) 



dVjn„,(f). V[r^+2n^(rO 



are toroidal and compressional operators, respectively. 
Note that both toroidal expressions ([^5)) and (P7)) involve 
the function r^~^^Y\f^{f) thus manifesting the relation be- 
tween TM and CM. In (1271) . the second term precisely 
gives the CM operator |31| . 

In Eq. p^ . the lowest-order fc-terms from ((23)) and 
([M)) exactly compensate each other and so we get 

Mvo. {EXp, k)=k [Mtor [EXp] -f A/com (^A/i)] . (29) 

Using the definition, Myor{EXp,k) = kM^or{EXp), we 
finally come to the relation 



M^oriEXp) = Mtor (EXp) + M^omiEXp) 



where 



(30) 



M^oriEXp) 
Mtor{EXp) 



M^om{EXp) 



2A + 1 



c(2A-h3) V A + 1 



d^nn«c(r)r^+inA+iM(70 



2c V 2A + 1 



d'rj„„,(r).r^+i 



1 A 



2c V A + 1 2A + 3 



dVr^+2l\;,^,(,^).(vxJ„„,(r)) 



A 



2c V 2A + 1 



d^r jnuc{r) ■ r 



A+1 



Yxx-l^,{r) - W ox^. ^ Yxx+it^i^ 



X 2X + 3 



-^^/^■V.-n,(.l (v. /.„.(.)) 



(31) 
(32) 



(33) 



Here, the TM and VM operators are the same as in comparison with the VM operator. Besides these forms 
but arc written in the forms convenient for the demonstrate the V x ]nuc{r) and V • ]nuc{r) origin of the 
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TM and CM operators, respectively. 

The expression for the VM operator ([51]) and the rela- 
tion (pO)) between VM, TM, and CM operators represent 
the main formal results of the present paper. 

Following (pO)) . the operators Mvor{EX^J,), Mtor{EXfi), 
and Mcom{EXiJ,) are of the same second order by k. They 
are given in (|31l) - (|33p in the current-dependent form. Us- 
ing the continuity equation (|2ip . the current-dependent 
CM operator ((33|) is straightforwardly transformed to the 
familiar density-dependent CM operator 



1 



2(2A + 3) 



d-'rpir^r^+'Yx^.ir) (34) 



as 



McomiEXfi) = -kACmiEXfi). 



(35) 



Note that the relation ([30)) requires the compensation 
of the terms ~ Y\A-iAt(r) the TM and CM operators. 
Thus a simultaneous use of these two operators is oblig- 
atory. The VM operator includes only Y\\+i^{'r) and so 
its matrix elements are determined by the current tran- 
sition density j\ a+i(''), as requested in It is easy to 
check that Myor{EXfi) gives the transition vorticity pT|) 
and so reproduces the results • 

The above formalism was derived for the case when 
the system is excited by the external electric field, i.e. 
for the excitation energy uj = hck > 0. The case of de- 
excitation is easily obtained by replacement /c — > — A: in 
the continuity equations ([5]) and (j2ip , k-dependent terms 
in (Uni), (HZl), dMI), and equations of Appendix C. The 
sign of Ms{EXijL, k) and density-dependent M{EX^) is 
changed as well. 



B. Dipole VM, TM, and CM operators 

The VM, TM, and CM are usually studied for the elec- 
tric /■^ = 1~ states Then the operators (|3ip - p3l) are 
reduced to 



MtoriEUi) = - 



2V3c 



d^rjnucir) 



(37) 



M^omiEllj) 



2^30 



d^rjnuc{r) 



(38) 



i'^r'Y^.^i^^ - (r^ - STMr\)Yio,i^)l 
[r'-dT,ol{r')or]Y^,{^^, 



(39) 



where (r^)Q = J d^rpo{r)r'^ / A is the ground-state 
squared radius. 

In (I37l)-(l39l), the terms Fiom(^> Yi^{?) include the 
center of mass corrections (c.m.c.) for T=0 excitations 
p^ . In the TM and current-dependent CM operators, 
the c.m.c. have the same magnitude. For the VM, the 
c.m.c. is zero, sec discussion in Appendix [BJ 

The expression in the square brackets of the TM oper- 
ator (1371) can be written as 



V X (fx V)(r3 - -(r2)„ T)Yi^{r^ , (40) 
which justifies a close relation between TM and CM. 



C. Discussion of VM, TM, and CM operators and 
vorticity criteria 

As shown above, the CM operator may be presented in 
the current-dependent psp and density-dependent (j34p 
forms. To the best of our knowledge, the former has not 
yet been used in the literature. Since the spin current jm 
is a curl of the magnetization, it does not contribute to 
the continuity equation and CM operator. Though the 
current-dependent form (j33p of CM formally involves j^, 

its contribution is annihilated by V • j-^uc (^) or the corre- 
sponding combinations of vector spherical harmonics. So 
actually both forms, ([33)) and ([34)) . of the CM operator 
do not depend on j„i • 

Both CM operators, ([33)) and ([34)) . are obtained 

through V • jnuc(f), which suggests their vorticity- free 
character. This is confirmed by the form of the CM ve- 
locity which, following the prescription [46l |. reads 



,(f)« Vr^+^yAM(f) 



(41) 



and so gives V x Wcom(f)=0. 

The current-dependent CM operator ([55)) includes the 

Y\\+ifj,{r) contribution which might be considered as an 
indicator of a vortical part. Indeed Y\\+ifj_{r) leads to 
the current component Jaa-i-i('') which, following Q, is 
responsible for the vorticity. However, for the CM this 
is misleading. The velocity ([^T)) can be straightforwardly 
cast into the form 



Vcor,i{r) OC 



A 



2A + 1 



(42) 



A + 1 2 ^ 
Yx\-i^,[r) - 2A + 3 "^^^+^^^^ 



whose curl is zero despite the Y\\j^i^{r) term. So. fol- 
lowing the HD criterion the appearance of a term 
Y\\+ifj,{r) in the mode operator and current is not yet a 
definitive signature of the vorticity. 

Altogether, we see apparent differences between two 
possible criteria for the vorticity: i) the HD condition 
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([2]) in terms of velocities, V x v{r) ^ 0, and ii) the condi- 
tion 01 in terms of transition current density components, 
3\ A+i('') 7^ 0. As shown above for CM example, a mode 
which is fully vorticity-free in the HD definition can have 
a substantial vorticity of the sort • The difference be- 
tween these criteria might be understood if we take into 
account that the vorticity density wwir) oc jAA-i-i('') 
was derived in Q, first of all, as a quantity completely 
unconstrained by the charge conservation rather than a 
purely vortical value in the HD sense. 

In fluid-dynamical and HD models [ill,[23, the velocity 
fields are chosen in form of the relevant external fields 
exciting the proper modes. Following this practice, the 
TM velocity reads 



vtorif) cx V X (fx \I)r^"''^Yxy.[r) 



2A + 1 



(A + l)(2A + 3)r 



A-l-l 



(43) 
(44) 



A 



A + 1 2A + 3 



>AA+1m(^ 



It is easy to check that V x wtor « r^Y\\^{r) and so TM 
carries the HD vorticity. 

The VM velocity constructed in the same manner is 



v^or{r) oc r^^^Yx-x^x^^{r) 



(45) 



It has the similar nonzero curl V r Ywuir). So 

the VM has the vorticity of both HD and Ref. Q sorts. 



IV. CALCULATION SCHEME 

The calculations of the excitation modes were per- 
formed within the separable random-phase approxima- 
tion (SRPA) model using the Skyrme energy functional 
[37I [38| . SRPA was earlier successfull y appli ed to descrip- 
tion of electric [37l - l4l| and magnetic |43l445j | GR in spher- 
ical and deformed nuclei. The approach was also used 
for the exploration of El strength near particle emission 
thresholds [42|. SRPA is fully self- consistent in the sense 
that both the static mean-field and factorized residual in- 
teraction are derived from the same Skyrme energy func- 
tional [49l - [5l| . The functional, £sk{p,T, J, j , s,T), in- 
cludes time-even (nucleon p, kinetic-energy r, spin-orbit 
J) and time-odd (current j, spin s, vector kinetic-energy 
T) densities. It also involves pairing (surface and vol- 
ume). Coulomb (direct and exchange), and c.m.c. terms 
[37I [ssl . [43j . The Galilean invariance of the functional is 
maintained in SRPA. The tensor spin-orbit contribution 
is involved through the squared spin-orbit densities ,7^. 
All the functional terms are kept in the mean field and 
residual interaction. 

The SRPA expands the RPA residual interaction self- 
consistently into a sum of separable terms, which dra- 
matically reduces the computational effort while keeping 
the accuracy of the full (non-separable) RPA [stI. l38j|. 



This makes SRPA extremely useful for systematic calcu- 
lations and tasks with very large configuration space, e.g. 
for description of GR in heavy deformed nuclei. 

For GR studies, the computational expense can be 
even more reduced by a direct evaluation of the strength 
function, thus avoiding the solution of RPA equations 
for a large multitude of individual states. The SRPA 
strength function works with the Lorentz weight and 
has very simple form [33, HI] . In the present study, the 
strength function for electric dipole modes reads 

Sa,iEl;uj)^i2nf ^ ^ |(«',|M„(^1m)|*o)PC(^^-c^.) 

(46) 

where 



C(w - w^) 



27r (^_^^)2 + 



(47) 



is the Lorentz weight with the smoothing width A 
and Ma{Elfi) is the electric dipole transition oper- 
ator whose type is determined by the index a = 
{vor, tor, com, com'}. Further, is the ground state, 

runs over the RPA spectrum with cigcn-frequcncies uj^ 
and eigen-states The Lorentz smoothing uses a 

width A=l MeV to simulate broadening effects beyond 
RPA (escape widths and coupling to complex configura- 
tions). The explicit form of the strength function ()46p . 
which does not directly involve the RPA solutions, is 
given elsewhere [13, HI] . 

In general, the VM, TM, and CM appear at many 
different multipolarities A^. Here we explore them in 
the 1^ = 1^ excitations of the doubly-magic spherical 
nucleus ^°^Pb. In this nucleus the pairing is absent. The 
calculations use the Skyrme parameterization SLy6 [s^ . 
which provides a satisfactory description of the E1(T=1) 
GR in heavy nuclei [4l|. For this parameterization, the 
tensor spin-orbit contribution is omitted. 

The calculations employ a cylindrical coordinate-space 
grid with the mesh size 0.7 fm. A large single-particle 
basis is used. The particle- hole 1~ pairs extend up to 
^ 65 MeV and, for E1(T=1) excitations, the energy- 
weighted sum rule with the isovector effective mass is 
exhausted by - 95% [4l|. 

The calculations involve both convection jc and mag- 
netization (spin) jm parts of the nuclear current, see 
Appendix [U] for more details. For jm, the isoscalar 
gl''P{T = 0) = {gP + gl')/2 = 0.88^ and isovector 
gl''P{T = 1) = (5? - gP)/2 = -4.70<; gyromagnetic 
factors are used, where gP = 5.58<j and 51" = — 3.82i; 
are bare proton and neutron g-factors and ^^=0.7 is a 
quenching parameter taking approximately into account 
the meson degrees of freedom [53| . It is easy to see that 
|.g^'P(r = 1)1 > \g^'P{T = 0)1 and so the main effect of 
the spin nuclear current is expected for T=l modes. 

The proton and neutron (q = n,p) RPA velocity fields 
arc determined in cylindrical coordinates {z,p) as 



Poiz,p) 



(48) 
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where 



208 



Jnuc 



) (49) 



is the current transition density for the RPA z^-state with 
the normalized forward and backward particle-hole (Iph) 



amplitudes 



and 



Further, p1^{z,p) is the pro- 



ton/neutron ground state density. 

The spurious cm. admixtures arc avoided by using the 
prescriptions from the Appendix|B] The SRPA equations 
and some important points, e.g. a choice of the generat- 
ing operators for the separable expansion, are sketched 
in the Appendix |DJ 



V. RESULTS AND DISCUSSION 

Results of the calculations for the nucleus ^^'^Pb are 
presented in Figs. [T]|9l 

In Figure E the VM, TM, and CM strengths in T=:0 
and 1 channels are compared. The strengths are com- 
puted for the transition operators (P5)) - (|55)) . For the CM, 
the current-dependent operator Mcom is used. Unlike its 
familiar density-dependent counterpart M[.^^, it has the 
same dimension as the VM and TM operators and so is 
more suitable for the comparison of the modes. All the 
modes are calculated with total nuclear current jnuc- 

Figure [1] shows that all the modes have basically two 
broad branches, a low-energy branch (LEB) at 5-20 MeV 
and a high-energy branch (HEB) at 25-40 MeV. The VM 
is well presented in both branches while TM and CM are 
mainly localized in LEB and HEB, respectively. Such a 
double-branch structure was found for the TM and CM 
in most of the previous theoretical studies 



15|. It is 

most probably related to El transitions with /S.N=1 and 
3 where N is the principle shell number. 

The double-branch structure of the El (T=0) strength 
was confirmed by various experiments [14l llSl . I20l425j . 
mainly in (a, a') scatte ring at small angles. The recent 
results of this reaction [24^[25| are depicted in Fig. [TJi). 
The reaction is considered as a common tool for measure- 
ments of the dipole CM(T=0) [H. 

As seen from the panel a) of Fig. [U our results do 
not reproduce the experimental energies and widths of 
the LEB and HEB. Neither of the mode centroids coin- 
cides with the experimental peak energies. Note that this 
is a common shortcoming of almost all theoretical stud- 
ies performed within various theoretical approaches (Tsj . 
Namely, the theory i) underestimates by 1-2 MeV the 
TM-dominated LEB and overestimates by about 4 MeV 
the CM-dominated HEB, ii) yields a much broader TM 
distribution and a too narrow CM one. The reason of 
the discrepancies is still unclear. Perhaps, this is partly 
caused by neglecting the coupling with complex configu- 
rations. Figure [1] shows that the VM and TM strengths 
are of the same order of magnitude in the LEB left flank 
for T=0 and in the whole LEB for T=:l (the difference 




10 20 30 

03 [MeV] 

FIG. 1: Isoscalar (T=0) and isovector (T=l) vortical, 
toroidal, and compression dipole modes in ^°*Pb, calculated 
with the SLy6 parameterization. The total nuclear current 
jnuc is used. The CM is computed with the operator Mcom 
from (|38[l . The hues with the arrows indicate widths and en- 
ergy centroids of the low-energy and high-energy branches of 
isoscalar El excitations observed in {a, a) reaction [2^.[25t. 



between T=0 and 1 cases is explained below in the dis- 
cussion of Figs. [2]and|3]). Perhaps, in these regions the 
TM is mainly vortical. We also see that VM dominates at 
the right LEB flank and is significant in HEB where the 
TM contribution is weak. The CM strictly dominates in 
HEB and has a noticeable tail in LEB at T=0. The latter 
is because of the coupling between TM and CM [H, Hlj . 
Obviously, the difference between VM, on the one hand, 
and TM and CM, on the other hand, is mainly explained 
by the terms ~ Yi^^ir) which are absent in Mvor but 
active in Mtor and Mcom- Just because of these terms. 
Fig. [1] cannot be used for a direct check of the relation 
([SO]) . For the same reason, the similarity of VM and CM 
strengths in the HEB cannot be considered as a signa- 
ture of the CM vorticity because VM and CM represent 
essentially different kinds of the motion, vortical versus 
irrotational. Instead, this rather means that both kinds 
of motion are presented by El AA''=3 transitions. 

The obtained results suggest that the VM(T=0) may 
be hopefully disentangled from other modes in (a, a') at 
the excitation energy ~ 16 MeV, where the VM(T=0) 
strictly dominates. For (e, e'), the pygmy region 7-10 
MeV seems to be most promising to observe VM and 
TM. In this region a dominant contribution ji2{'r) of the 
nuclear current is expected. 

In Figure [21 the isoscalar VM and TM strengths, cal- 
culated with the complete jnuc = jc + jm, convection jc, 
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FIG. 2: Calculated isoscalar (T=0) vortical, toroidal, and 
compression dipole modes in ^"^^Pb. For the VM and TM, 
the strengths with the total jnuc, convection jc, and magne- 
tization jm current contributions to the transition operators 
are shown. 



and magnetization j,„ nuclear currents, are compared. It 
is seen that the jm contribution is weak and so the T=0 
VM and TM are mainly of the convection nature. This 
is especially the case for the HEB where the jm contri- 
bution, being mainly of low-energy spin-flip character, 
vanishes at all. The weakness of the jm weight in T=0 
channel is naturally explained by the low values of the 
gyromagnetic factors g](T = 0), as mentioned in Sec. 
TV\ Fig. [5] also exhibits the CM. Unlike Fig. [H here 
the familiar density-dependent operator M^^„^ is used. 
Following (pS)) . Mcom is less energy- weighted than Mcorn 
and so gives a more comparable CM strength in LEB 
and HEB. As discussed above, the CM is determined by 
V • jnuc and so is purely irrotational. It has no any con- 
tribution from jm and thus is fully of convective. 

In Figure O the VM, TM, and CM are shown in the 
T=l channel. As compared to the previous T=0 case, 
we see dramatic changes in the magnitude and composi- 
tion of VM and TM. In the LEB, these modes become 
stronger and dominated by the jm contribution. The rea- 
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FIG. 3: The same as in Fig. [2] but for the T=l modes. 



son of the changes is obvious. The isovector spin factors, 
g"'^{T = 1) = — 4.70i;, are much larger than the isoscalar 
ones, 5^P(r = 0) = 0.88<;. So, the T==l spin contribu- 
tion grows about {g^'P{T = l)/gl''P{T = 0))^ ~ 29 times. 
It becomes dominant and significantly increases the to- 
tal VM and TM strengths. Note that this effect does 
not concern the HEB which remains purely convective. 
Besides, the jm effect is zero for the CM. 

The next point to be considered is collectivity of the 
modes. To this end. Figs. S] and [S] compare RPA and un- 
perturbed particle-hole (Iph) strengths. It is seen that 
the RPA residual interaction noticeably down-shifts the 
strength for T=0 and up-shifts it for T=l. The maximal 
collective shifts (defined as a difference between RPA and 
Iph peaks) take place in the CM, where they reach 1-2 
MeV in LEB and 2-4 MeV in HEB. The HEB shift is com- 
parable to that of the E1(T=1) giant dipole resonance 
(GDR) exhibited in Fig. [5] and so is indeed very large. 
This indicates that HEB modes, VM and CM, are collec- 
tive. The LEB modes, for exception of a few high peaks, 
are less collective. The LEB almost coincides with the 
region of the unperturbed Iph dipole strength depicted 
in Fig. [H] and so for the LEB the single-particle aspect is 
also important. These observations are confirmed by in- 
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FIG. 4: Comparison of SRPA and Iph strengths for the T=0 
modes. For VM and TM, the strengths are computed with 
the total nuclear current jnuc- 



spection of the detailed structure of the RPA states and 
agree with the previous studies [l^ for the high-energy 
CM and low-energy TM. 

The next figures take a closer view at the detailed 
structures of a few most collective LEB modes. In Figs. 
ElEl the VM, TM, and CM neutron and proton veloc- 
ity fields dHj) for particular T=0 and T=l RPA = 0) 
states are shown. The states in the region 6-10 MeV with 
the maximal VM, TM, and CM responses are considered: 
8.3 McV (VM), 8.7 McV (TM), 7.1 McV (CM) for T=0 
and 9.1 McV (VM), 9.8 McV (TM), 8.8 McV (CM) for 
T=l. These states are easily recognized in Figs. |2]|S]as 
highest peaks at the left side of the LEB. The states com- 
bine collectivity and single-particle effects: their struc- 
ture is a coherent superposition of many Iph components 
with maximal contributions 20-40%. The velocities are 
plotted in cylindrical coordinates (z,p) and, for simplic- 
ity, only the quadrant (z > 0, p > 0) is presented. 

Figures [7][9| show that velocity fields are rather in- 
volved, indicating a complex structure of the RPA states. 
The clear imprints of the vortex motion are seen as lo- 
cal curls. Sometimes, the fields well reproduce a typical 
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FIG. 5: The same as in Fig. |4]but for the T=l modes. 
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FIG. 6: SRPA and Iph strengths for the E1(T=1) GDR. Like 
for other modes, the strength functions are plotted without 
the energy weight. The experimental width and energy [sj] 
are shown by the horizontal line and arrow, respectively. 



toroidal picture, see e.g. the T=0 proton velocities for 
VM and CM in panels b) of Figs. [TjandEl The curls for 
the CM may be explained by its strong coupling to TM 
in LEB. In Fig. |9^), the strong dipole component is also 
seen. However, in most of the panels, a large impact of 
the single-particle motion, e.g. of the nodal structure of 
the leading Iph components, distorts the familiar coUcc- 
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FIG. 7: a)-b): Neutron and proton vortical velocity fields 
v^{z,p) for the state oJv ~ 8.3 MeV with a maximal T=0 
vortical response. c)-d): The same for the state uji, = 9.1 
MeV with a maximal T=l vortical response. For a better 
view, the velocities are amplified by the factors 50 (a-b) and 
20 (c-d). 



FIG. 9: The same as in Fig. [7] but for the states uu^ = 
7.1 MeV (T=0) and uj^ = 8.8 MeV (T=l) with a maximal 
compression strength in the low-energy branch. 



VI. CONCLUSIONS 




FIG. 8: The same as in Fig. [7] but for the states lui, — 
8.7 MeV (T=0) and uj^ = 9.8 MeV (T=l) with a maximal 
toroidal strength. 



tive TM and CM flows and considerably complicates the 
picture. 



The nuclear vorticity and relevant vortical, toroidal, 
and compression modes (VM, TM, and CM) were ex- 
plored on general formal grounds. The operators of all 
three modes were derived as second-order terms in the 
long- wavelength expansion of the electrical multipole op- 
erator and its subsidiary counterpart built following the 
concept 0. For the first time, the vortical operator was 
constructed and related to its toroidal and compression 
counterparts. The detailed comparison of the modes and 
their operators was done. It was explicitly shown that, 
while VM (CM) is vortical (irrotational) by construction, 
the TM is of a mixed character. The vorticity criteria 
from HD and arguments based on the decoupling to the 
charge conservation Q were inspected. It was shown that 
the latter deviates from the HD definition and may lead 
occasionally to misleading conclusions as, e.g., a vorticity 
of the CM. 

The electric dipole VM, TM, and CM were com- 
puted and investigated in ^°®Pb within the self-consistent 
Skyrme random-phase-approximation (RPA) approach. 
Unlike most of the previous studies, both convection 
and magnetization (spin) parts of the nuclear current 
were taken into account and both isoscalar (T=0) and 
isovector (T=l) channels of the modes were analyzed. 
It was shown that VM and CM have low-energy and 
high-energy branches while TM mainly appears in the 
low-energy branch. The CM strictly dominates in a high 
energy branch. In the T=0 channel, the VM and TM 
are almost completely determined by the nuclear convec- 
tion current while in the T=l channel, their low-energy 
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branches are strictly dominated by the spin current. This 
strong isospin effect is straightforwardly explained by low 
(high) values of T=:0 (T=l) spin g-factors, which drasti- 
cally changes the ratio between the convective and spin 
contributions of the current. The effect cannot appear in 
the CM since this irrotational mode has vanishing spin 
contribution. 

The collectivity was found strong for the high-energy 
VM and CM and rather weak for the low-energy VM, 
TM, and CM. In the latter case, the velocity fields of the 
modes are rather involved. This is partly a consequence 
of the complex structure of the RPA states mixing the 
collective with detailed Iph contributions. This holds in 
particular for the vorticity which seems to be of both 
single-particle and collective origin. 

The VM, TM, and CM were shown to be closely related 
and, at the same time, display considerable differences in 
their detailed strength distribution. These modes seem 
to represent one family with complementing aspects. It 
would be interesting to analyze the results of our study 
more deeply and use them to disentangle the modes in 
the (e, e') and hadron reactions. This is in our next plans. 
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Appendix A: Removal of the charge conservation 
constraint 

In Sec. Ill A, the vortical operator M^or{E\^, fc), com- 
pletely unconstrained by the continuity equation, is built 
from the M{E\ijl, k) operator by replacing the curl of the 

nuclear current V x inuc{f) by the vortical quantity 

t^A = [V X jnuc{r)\ \ - jkc [V/5(r) x r\x (Al) 

where the terms [...]a assume A-components of the mul- 
tipole expansion of the values inside the brackets. Below 
we present arguments which motivate (jAip and compare 
it with the prescription Q. For the simplicity, we con- 
sider only the convection current and neglect isospin. 



1. Simple arguments 

Using the HD definition of the velocity field ([3]), we 
may write the truly vortical quantity 

Po{r)y xv^{f) = V X 5>(r) - V po{r} x il^{r) . (A2) 

Except for the second velocity-dependent term, the r.h.s. 
of (|A2|1 reminds the operator construction (|Aip and thus 
may be used for justification of uj\ as a vortical quantity. 

To make (|A2p closer to (|A1I) . it is worth to express r 
in terms of the (global) velocity operator 

V = r = —[H ^ ikcf . (A3) 

This suggests the replacement 

V — )■ ikcf (A-4) 

in (IA2[) and thus casts (jA2p to a form similar to (jAl[) (up 
to the multiplier 1/A). Note, however, that the velocities 
in (jA2p and (|A3[) are not the same. The replacement 
of the velocity field v^i'^) by a global velocity becomes 
strictly valid only in the sum-rule limit when all multi- 
pole strength is concentrated in one collective state. It 
remains probably an acceptable step for distributed spec- 
tra which often gather around a strongly collective mode. 

2. Correspondence of the recipes 

It is also worth to relate the recipe (|T0)) for the transi- 
tion densities (3] with our recipe ([T5|) and its analog (jAl|) 
for the operators. As was mentioned in Sec. Ill A, both 
recipes serve to build the vortical quantities, though by 
different ways. While ([TO)) excludes from the current the 
jxx-i{r) terms, the recipe (|15p leads to exact compensa- 
tion of the lowest-order fc-terms in the vortical operator. 

One may show that the recipes ((TU)) and actually 
correspond each other. This may be done by treating 
(|15p in terms of the transition densities © and currents 
(O. Using the relation between the density and current 
expansion multipoles 

kcpxir)^ -/^(|:-^)..AA-.(0 (A5) 
/ A+ 1 , d A+2_ , , 

one may show that (fTS]) indeed leads to the unconstrained 
vortical transition density (|10p . In this case, the second 
r.h.s. term of Eq. (|Aip has the form 

zkc [V%.(r) X r] = ^ a{;5AA(r) Yxx^{^ (A6) 

Xfj, 

with the miihipoles 

Sxx{r) = ^X{X + l)kcpx{r) . (A7) 

Being scaled by 1/A, these multipoles coincide with the 
expansion coefficients of 6Sfi{r) in 
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Appendix B: Extraction of spurious admixtures 

The isoscalar electric dipole VM, TM, and CM may 
have spurious admixtures related to the center of mass 
motion of the nucleus. There are various methods to 
derive the center of mass corrections (c.m.c), see e.g. 
[1, [m, [ssl - fSTi . Most of the methods neglect the velocity- 
, spin-, and spin-orbit-dcpcndcnt terms in the nuclear 
interaction and assume a closure relation where all exci- 
tation is contained in one single collective state. Hence 
the methods are approximate. Nevertheless, being sim- 
ple and convenient, they are widely used in the Skyrme- 
Hartree-Fock calculations [iBl . Below we use the method 
[56j to derive the c.m.c. for isoscalar TM and CM. For 
the isoscalar VM, the c.m.c. is shown to be zero. 

Let's consider a system with the Hamiltonian H = 
T + V whose interaction V does not depend on spin and 
velocities. Then, in the isoscalar (T=:0) case, for any 
one-body external field F = firl), where /(n) is an 
arbitrary function of nuclcon coordinates, one may write 
the sum rules WOl 



Y,ujAO\p{r)W){^\F\0) 

V 

j2{o\j{r)W){'^\m 



1 

1 

2mi 



V • [po(rlV/(rO](pi) 
Po(r)V/(f) (B2) 



for the isoscalar transition density (0|/5(r)|:/) and current 

{0\j{r)\v). The sum runs through the full set of the exci- 
tation eigenstates |^) {H\i^) = uj,y\i/)) with the eigenval- 
ues w^. Further, poir) = (0| S{f—fi)\0) is the nuclear 
ground state density, and m is the nucleon mass. 

If the external field mainly excites one single collective 
state P, then only the term (i/ = P) survives in (jBip - (jB2p 
and the corresponding transition density and current are 
uniquely specified through po{r) and /(f) (46j : 



{0\p{r}\i? 



1 1 



2ma;p(zy|i^|0) 



2mi {y\F\G) 



(B4) 



Hence we get simple recipes for the transition densities 
and currents, 



(B5) 
(B6) 



8p[r) cx V • [po(r1V/(rl] , 
6j(r) (X po{r)'^.f{'r) cx po{r)v{r) 



to be used in further c.m.c. derivation. Here (jB6[) 
actually defines an irrotational flow with the velocity 
v{r) ^ V/(f). 

A change of the expectation value (O) of any one-body 



operator O = J2i 0(^1); caused by an external field F, is 
S{d) = J d^r 5p[r) o(f) (B7) 
cx j d^r o(r) V • [po(f) V/(f)] = (B8) 
= - I d\ po(f) V/(f) • Vo(f) (B9) 
cx - / d^r Sj{r} ■ Vo(f) (BIO) 



where we use the relations (|B5[) - (jB6[) . Both (5p(f) 
and Sjif) may be applied to determine 5{0). These 
cases are suitable for the modes determined by the 
density-dependent and current-dependent operators, re- 
spectively. Note that Eq. (|B10p with Sj{f) is general 
and may be derived and used by itself, rcgardlessly to 
the formalism (|Bip - (jB9p and character of the flow. The 
quantity Sj^r) is then determined by Eq. ([3]) and, un- 
like (jB6p . the velocity v{f) of the flow can be not only 
irrotational but also vortical or mixed. Altogether, the 
relations (jB7[) - (|B10[) may be reduced to an expression 



HO) 



d^r po{r) v{r) ■ Vo(r') 



(Bll) 



with v(r) covering both irrotational (|B6I) and more gen- 
eral cases. 

If (O) is invariant with respect to the field F, then the 
requirement 5{0) = is kept and (jBlip gives 



For 



d^r po{f) v{f) ■ Vo(f) = 



o{r) = rY*ir), 



(B12) 



(B13) 



the operator O = J^i o(fi) describes the cm. coordinate 
and (|B12p means that this coordinate is not affected by 
the field F. By using 

Vo(f) = VrY{^i^ = VS Y.l^i?}, (B14) 

the conditions (|B12p is cast into 



d^rpoif}vir}-Y*,^ir) =0 



(B15) 



The condition (jBlSp may be applied to the density- 
dependent CM operator M^^j^{Elp) by imposing the ir- 
rotational velocity 

Vco,n' (r) = V/co,„ (f) = VYi^ (rO (r^ - rjr) (B16) 

where the second term is the actual c.m.c. with the co- 
efficient 77 to be determined from (|B15p . Then, assuming 
spherical nuclei (po(f) = Poi'i')) and using 



(B17) 
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we get 



with 



»7 3 ;o 



t"^ Pf){r)dr 



(B18) 



(B19) 



and finally the expression (P^)) for M^^j^{Elfi). 

For the current-dependent VM, TM, and CM operators 
dMll-dSH]), we put to (|B15P the velocities 



Vvorir) = 




- ■qYwt.ir) , 




(B20) 


Vtorir) = 


—r Y,,, 


{^+Y,^^.{^{r^ 




, (B21) 


^com(^) — 


^/2 2v 


.i^-Y,o,{f)ir' 




, (B22) 



where the second terms with 77 are the relevant c.m.c. 
For exception of Vcom i'^) , these velocities are not reduced 
to the gradient form. Actually, they are taken in the form 
of the external fields involved in the operators p6p - ([55)) 
and generating the corresponding modes. Such presen- 
tation is in accordance with the self-consistent treatment 
of nuclear excitations [l[, which is done here in terms of 
small variations Sj{r) of the nuclear current. 

By using (jB20p - (|B22|) . the requirement (jBlsp gives 



r] = 



(B23) 



Appendix C: Nuclear density and current operators 

The density operator reads 

E e%J2(^ir-rk) (CI) 

q=n,p keq 

where e^g are proton and neutron effective charges. 

The operator of the full nuclear current consists of the 
convective and magnetic (spin) parts [ist 



jnucir) = jcir)+j,n{r) 
where 



eh 



q—n^p 



keq 



frr 



(^) = y E"^ ^ Kk5{r-rk), 

keq 



(C4) 



and Sq is the spin operator, /xjv is the nuclear magneton, 
g1 is the spin g-factor, k numerates the nucleons. 
The T=0 modes use the values 

1, 5^'>f(r = 0) = i(g^'+5?) (C5) 



while the T=l modes employ 



-1, g-^'P{T=l)^-{g:-gl). (C6) 



for VM and 

{r')o (B24) 

for TM and CM. Thus we get the corrected expression 
(1271)- (123) for Mtor{Elfi) and Mcom{Elfi). Note that vec- 
tor harmonics Yi2fi{r) related to the vorticity do not con- 
tribute to the c.m.c. This reflects the physical fact that 
vorticity, being a curl flow, must be fully decoupled from 
the cm. translation motion. Hence the c.m.c. is zero for 
the VM. On the other hand, the TM is not completely 
vortical and so its c.m.c. does not vanish. 

Note that the above c.m.c. are approximate. Indeed, 
the calculations give for the VM, TM, and CM responses 
two broad structures, which actually do not meet the 
sum-rule condition of excitation of a single collective 
state. Besides, the prescription (46j uses the commutator 
[H, p] where H is assumed not to include the terms with 
velocity-, spin-, and spin-orbit dependence. However, the 
effect of spin-dependent terms in the commutator is ob- 
viously zero. The momentum (velocity)-dependent inter- 
action does not matter for the Galilean-invariant Skyrme 
functional (the most common case) but may be impor- 
tant if this invariance is violated. The spin-orbit interac- 
tion may affect the c.m.c. 



Appendix D: SRPA equations and generator 
operators 

The SRPA Hamiltonian is self-consistently derived [13, 
Issl l from the functional 

£ = £kin + £sk + ^^pair + ^^Coul (Dl) 

involving kinetic-energy, Skyrme, pairing and Coulomb 
terms. The Skyrme functional fsk(p, t, J, j, s, T) de- 
pends on time-even (nucleon p, kinetic-energy r, spin- 
orbit J) and time-odd (current j, spin s, vector kinetic- 
energy T) densities. The Hamiltonian reads [13, HI] 

H - ZiHFB + Kcs (D2) 

where Ahfb is the HFB mean field 

and Vres is the separable residual interaction 
1 ^ 

^'^'''=2 ^ i'^kk'XkXk' +Vkk'YkYk') (D4) 

k,k' = l 
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with one-body operators 
Yk = ij j d'rd'r' sjfsj . ^tQ'-"^"-])'^"'^ ' 

a_,a'_ ~ ~ 

and inverse strength matrices 

Kll,^^i{[Pu,Xk>]) , Vkk' ^ ~^{[Qk,%]) ■ (D5) 

Here a+and a_ enumerate time-even J^^ and time-odd 
Jq_ densities, respectively; Jq,^ are the density operators; 
Qk and Pk ~ i[H, Qk] are time-even and time-odd hermi- 
tian generator operators. The operators of the residual 
interaction Xj- and Yk are time-even and time-odd by 
construction, respectively. 

The single-particle Hamiltonian /ihfb is determined by 
the first functional derivatives of the initial functional 
(|Dip while operators Xk and Yk are driven by the sec- 
ond functional derivatives of the same functional. The 
residual interaction includes all the possible terms aris- 
ing from (|D1[) . Hence the model is fully self-consistent. 
The number K of separable terms in (|D4I) is determined 
by the number of the generator (input) operators Qk- 
Usually we have K = 3 — 5. This results in a low rank of 
the RPA matrix and so in an efficient calculation scheme. 

The SRPA formalism itself does not prescribe the form 
of the generators Qk and Pk- At the same time, their 
choice is important for a fast converge of the factorized 
residual interaction V^cs to the true one with a minimal 
number of separable terms. The set of the generating 
operators is introduced so as to initiate in the nucleus all 
the relevenat kinds of motion for the considered modes. 
For time-even modes, the initial generators Qk are chosen 
first and then their time-odd counterparts are determined 
from Pk ~ i[P[ , Qk]- Instead for time-odd modes, the ini- 
tial generators Pk arc inserted and then their time-even 
counterparts Q^ — i[H,Pk] are determined. The genera- 
tors may be arbitrarily and separately scaled, which docs 



not influence the results. The coupling of the modes (e.g. 
of electric and magnetic ones in deformed nuclei) may re- 
quest both time-even and time-odd generators in the set. 
The optimal sets of the generators were developed for 
E1(T=1) [13, Hi] and spin-flip Ml GR. 

Here we use the minimal sets of the generators suitable 
for the description of VM, TM, and CM. The generators 
cover the main parts of the corresponding operators and 
take into account the coupling between the modes. For 
the VM, they are 

A = y" dVje(r)r2yi2^(f), 
A = y d^rl^ifyY,2^.ir^, 

Q3 = Jd'rpir)[r^~^{r^),r]Y,^{r), (D6) 

i.e. cover the time-odd parts of the vortical operator 
with the convection and magnetization currents as well 
as the time-even compression operator with the c.m.c. 
(to prevent generation of the spurious motion). 

For the TM, the generators read 

A = 1 d'r lir)[r^^Y^,,i^^ + Y,,,{?){r^ - (r^)^], 
A = 1 d\ J^(f)[r2^fi2^(r^ + yiop(r1(r2 - (r^)^], 
Q3 = 1 d\p{r)[r^ - ^(r2)„r]ri„(r^, (D7) 

i.e. cover the time-odd parts of the toroidal operator 
with the convection and magnetization currents as well 
as the time-even compression operator. The TM needs 
the c.m.c. and so now this correction is included to all 
the generators. 

Finally the set for the CM includes the compression 
operator itself and the convective toroidal generator 

A- The generators with the magnetic current Jm(^) are 
not involved since their effect on the CM is zero. 
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